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Abstract. 

Loss of atoms from a *^Rb condensate on passage through a Feshbach 
resonance is analyzed using the generalized parametric approximation that takes 
into account quantum many-body effects. These effects lead to a substantial 
increase of the losses. A better agreement with experiments is achieved, compared 
to predictions of mean-field theories. The method provides much insight into the 
quantum effects involved, and on the nature of entangled atom pairs produced by 
the loss. 
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The presence of a Feshbach resonance is supposed to provide a tool for controlling 
the interatomic interactions determining the properties of a Bose-Einstein condensate 
(BEC) (see ^ and references therein). A Feshbach resonance occurs when the energy 
of a pair of atoms in the condensate is close to that of a metastable molecular state. 
The scattering length, as a measure of the mean interatomic interaction, varies strongly 
as a function of the energy mismatch between the two states. This mismatch can be 
controlled by applying a varying magnetic field. The energies of the two states can be 
brought closer to each other, as the two states have different Zeeman shifts. 

The effect was studied first in the MIT experiment on Na by applying a 
time- varying magnetic field B = B (t). The experiment resulted in a large condensate 
population loss. In order to provide a theoretical explanation of these experimental 
results, two mechanisms have been suggested. The first one is a collisional deactivation 
process relating the loss to atom-molecule and molecule- molecule inelastic 
collisions. This mechanism has been analyzed in ^ |3] using coupled mean- field 
(MF) equations for the atomic and molecular condensates. The second mechanism 
is an excitation process involving a crossing of the resonant molecular state into 
non-condensate atomic states. It has been analyzed at first in 0] as a dissociation 
of single molecules, without taking into account many-body effects. The combined 
effect of both mechanisms was studied in [5] , where the crossing mechanism has been 
incorporated into MF equations by introducing a width to the molecular condensate 
state. It has been shown that in the case of the MIT experiment both mechanisms 
contribute to the loss comparably and non-additively. The crossing loss mechanism 
and non-condensate states are especially important in the case of the extremely strong 
Feshbach resonance in ^^Rb studied in the more recent JILA experiments 0. 
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The non-condensate atoms are formed as entangled pairs in two-mode squeezed 
states [7]. Several methods have been suggested earlier, allowing a correct treatment 
of such essentially quantum states. A numerical solution of stochastic differential 
equations in the positive-P representation has been used for this purpose in 
In the Hartree-Fock-Bogoliubov formalism (HFB) |^ ^] the MF equations were 
complemented by equations for the normal and the anomalous densities, describing the 
second-order correlations of the non-condensate atomic fields. These correlations are 
taken into account to the same accuracy in the parametric approximation [3 1111 112| 
and in the microscopic quantum dynamics approach |13[I14| . Some qualitative results 
have been also presented in ^^l- A modified version of the MF theory ^Bj differs 
from the HFB method by neglecting the normal density, which actually reduces the 
problem to a two-body one (see ^3E|)j removing quantum many-body effects. An 
attempt to describe the JILA *^Rb experiments using this theory in has therefore 
met with little success. 

The present work uses the parametric approximation T in order to describe the 
JILA experiments jB| . The incorporation of the quantum effects leads to a much better 
agreement with the experimental results compared to MF theories. The results have 
been preliminarily reported in |19| . A more recent work j2()j , performed independently, 
uses the microscopic quantum- field-theoretic approach of jl4| . in order to obtain 
the same goal. The approach of [201 takes account of the spatial inhomogeneity of 
the condensate, and therefore yields a better agreement with the experimental data, 
compared to the one presented here. 

The approach of jj^ ^sed by [201 , based on the numerical solution of a nonlinear 
partial-integro-differential equation, may provide the ultimate word concerning 
accuracy. But it lacks the transparency of the parametric approximation as used 
here. At the cost of avoiding inhomogeneities, the present method provides a greater 
qualitative insight, such as the occupation of the various non-condensate states, from 
which useful information can be derived regarding the profile of the entangled atom 
pair production. Also one can easily trace causal effects, such as the relation between 
the the excessive condensate loss and the effects of quantum Bose enhancement on the 
curve-crossing process (see QH). 

Following the generalized parametric approximation |7] , let us consider a system 
of coupled atomic and molecular fields described by annihilation operators in the 
momentum representation (p,*) and {p,t), respectively. The coupling of the 
atomic and molecular fields contains a product of two atomic creation operators and 
therefore describes the formation of entangled atomic pairs. Spatial inhomogeneity 
due to the trapping potential and the effects of elastic collisions are neglected here. 

Let the initial state of the atomic field a,t t — to he a coherent state of zero kinetic 
energy 

(P, to) |in) - (27r)3/2 ip, (to) S (p) |in}, (1) 

where |(pa (to) P — "-o (to) is the initial atomic condensate density and |in) is the 
time-independent state vector in the Heisenberg representation. A pair of condensate 
atoms forms a molecule of zero kinetic energy. Therefore the resonant molecules can 
be represented by a mean field ipm (t) as 

(in|*,„ (p, t) |in) = (27r)3/2 ^ ^ (3) 

where \(pm (t) P = rim (t) is the molecular condensate density. Fluctuations of the 
molecular field due to collisions involving non-condensate atoms are neglected. 
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The outcome of atom-molecule and molecule-molecule deactivating collisions is 
introduced by adding molecular "dump" states. The elimination of these states in a 
second-quantized description leads to the equation of motion for the atomic field J, 

„2 



iP,t) = 



\-<^a (t) - h\(Pm (t) I 

Zm 



(P, i) + 2.9>„ (t) (-P, t)+iF (p, t) , (3) 



where m is the atomic mass, ea{t) — — (^) ^ Bq) is the time-dependent Zeeman 

shift of the atom in an external magnetic field B (t) , relative to half the energy 
of the molecular state, ii is the difference in magnetic momenta of an atomic pair 
and a molecule, and Bq is the resonance value of B. The coupling g of the atomic 
and the molecular fields is related to the phenomenological resonance strength A as 
\g\'^ = 27r|aa/^|A/TO jSI, where aa is the background elastic scattering length. The 
quantum noise source F (p, t) and the parameter 7 describe the effect of deactivating 
collisions. The deactivation plays important role in general case and is included 
in the present calculations. However for the conditions of experiments [S] the results 
of calculations are insensitive to these processes (a variation of the deactivation rate 
coefficients from to 10~^ cm'^/s change the results by less than 1%). The contribution 
of deactivating collisions is therefore neglected in following analysis. The atomic field 
operator can then be written in the form 

(p, t) = (p, to) (p, t) + *1 (p, to) V. (P, t) . (4) 

The c-number functions ^Jc,s {p, t) are solutions of the equations 



iij^cs (p, t) 



r „2 



*,.(?,<) + 29 V™ (')viap.i). (5) 



given the initial conditions ipc {p, to) = li V's (p, ^o) = 0, and the constraint 

|^aP,i)l'-|V^.(p,i)l' = l- (6) 
The two-atom correlation functions are expressed in terms of these solutions as 

(inl^t (p, t) (p', t) |in) = (27r)^ no (t) 5 (p) 5 (p') + n. (p, t) (5 (p - p') (7) 

(inl^-a (p, t) ^-a (p', t) |in) = {2T:f mo (t) 5 (p) 5 (p') + m, (p, t) 5 (p + p') , (8) 
where 

no (t) = 1^, (t) p (9) 
is the condensate density, 

(t) = (in|*a (0, t) |in) = Vc (0, t) (^^ (to) + V-s (0, t) (to) (10) 
is the atomic mean field, 

n,(p,t) = |V'.(p,t)|2 (11) 
is the momentum spectrum of the non-condensate atoms, and 

mo (t) = (t) , (p, t) = V. {p. t) V-c (p, t) (12) 

are the anomalous densities of the condensate and non-condensate atoms, respectively. 
The atomic density 

Ua (t) = (27r)~' j d^pid^P2 exp [i (p2 - Pi) • r] (in|*t (p^, t) #a (P2, t) |in), (13) 

X A system of units in which Planck's constant is ^ = 1 is used below. 
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then appears to be r- independent, and comprises the sum {t) — uq {t)+ J (PpUg (p, t) 
of the densities of condensate atoms and of non-condensate (entangled) atoms in a 
wide spectrum of kinetic energies. 

The equation of motion for the molecular mean field ipm (t), obtained by the 
removal of the effects of deactivating collisions, followed by a mean-field averaging, is 

_9_ 

27r2 



i(Pm (t) = gmo (t) + / p^dpms {p, t) , (14) 



where the cutoff Pc in momentum space is necessary in order to avoid divergences. 
The final results, after a renormalization of the detuning ea (t) in a manner similar to 
the one used in 9 , are insensitive to the value of the cutoff momentum if it is large 
enough, so that p^/ (2m) > max(— {t)). 

Equations Q, 10), 110() . 1)111) . and (|12|l lead to the following equations of motion 
for the atomic mean field, and the normal and anomalous densities of the non- 
condensate atoms, 

i^a (t) = €a (t) ipa (t) + 2fJ*ipl {t) ifim (t) (15) 

ifis {p, t) = 2g>™ (t) m* {p, t) - 2gLp*^ (t) m,, {p, t) (16) 



inis {p, t) 



- + 2ea (t) 



m, {p, t) + 25 V™ (t) [1 + 2n, (p, t)] . (17) 



These equations, in combination with (|14|l are the momentum representation of the 
HFB equations used in Therefore the present approach becomes mathematically 
equivalent to the HFB one having, however, some advantages in simplifying the 
numerical calculations. 

The results presented here were obtained by a numerical solution of Eqs. Q on 
a grid of values of p, combined with Eq. ((Tljl . The values of A = 11 G, Bq — 154.9 G, 
and loal = 450 (in atomic units) are taken from [2], and the value oi ^ — —2.23 (in 
Bohr magnetons) is taken from ^U] . The dynamics is calculated for a linear variation 
of the magnetic field B = Bq — Bt from 162.3 G to 132 G with various sweep rates B, 
chosen in accordance with the experiments 

The results are sensitive to the initial value of the magnetic field. This sensitivity 
can be attributed to the large crossing width Tcr which substantially exceeds the 
range of the detuning ea {t) along the full sweep. Therefore, in some sense, the system 
is always within the resonance. As the initial values of (/?„i, ipci and ips were substituted 
the steady-state values calculated using the initial magnetic field as a constant. The 
calculations were performed for a homogeneous BEC with the initial density 10^^ 
cm~^. 

An example of the resulting dynamics is presented in figure ^ It demonstrates 
that the main mechanism of the condensate loss is a conversion to non-condensate 
atoms by the crossing mechanism, while the molecular occupation, and hence the 
deactivation losses, are negligibly small. The effect of Bose enhancement on the 
crossing from the molecular state to non-condensate modes is proportional to the 
mode occupation presented in figure QJj. 

Figure EJl compares the results of calculations to the experimental data 0. 
The results of mean-field calculations using the method are presented as well. 
It demonstrates the better agreement with the experimental data attained by the 
parametric approximation, compared to the MF results. The drastic increase of 
loss at the slower sweeps is actually related to the effect of Bose enhancement, as 
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Figure 1. (a) Time dependence of the densities of the atomic condensate ( ), 

the molecular condensate ( ), non-condensate atoms ( — • — ) and of the 

mean squeezing parameter f I18i ( ), calculated for the ramp speed B = 5 

G/ms. (b) The occupation ris {p,t) of several non-condensate modes with energies 
p2/ (2m) = 1.1 nK ( ), 4.4 nK ( ), 9.8 nK { ), and 17.5 nK ( ). 




Figure 2. (a) Ratio of the surviving atomic condensate density no (t) to the 

initial one no (to) ( ), and the maximal non-condensate state occupation 

( ), as functions of the inverse ramp speed, calculated using the generalized 

parametric approximation, in comparison with the experimental data by Cornish 
et al. E] (• ), and the results of mean-field calculations based on (5] ( — ■ — ). (b) 
Energy spectra of non-condensate atoms Us {E, t) I19i ( — ■ — ) and the squeezing 

parameter r{E,t) 1181 ( ), calculated for the ramp speed B = 5 G/ms at the 

time t = 1 ms. 



demonstrated by the plot for the maximal value of the non-condensate state occupation 
ris in figure The disagreement between the parametric approximation and the 
experimental data at intermediate sweep rates is related to the effects of spatial 
inhomogeneity neglected in both theories, but taken account of in the calculations 
120. 

The non-condensate atoms are formed in two-mode squeezed states, which 
are entangled [7]- Entangled atoms can have many applications in quantum 
measurements, calculations, and communications (see and references therein). 
The amount of entanglement can be measured by the energy-dependent squeezing 
parameter r{E,t) and a mean squeezing parameter f{t) [7|, 

.J, ,^ 1, l + 2n,{p,t) + 2\fnAp,t)\ J d^pnsip,t)r{£^,t) 

'^^''''^ A'''\l + 2n,ip,t)-2\ms{p,t)\\' "^^'^ J d^pn^ip^t) ^'^^ 

where E = p^ j (2m) is a kinetic energy of non-condensate atoms. In our calculations 
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the maximal squeezing takes place at the slowest sweep. Even in this case the 
mean squeezing does not exceed the value of 1.1 (see figure^), and the energy- 
dependent squeezing does not exceed the value of 1.4 (see figurel^b). These values are 
substantially less than maxf(i) « 3.1 and ma.xr{E,t) « 3.6 calculated for the weak 
resonance (A w 9.5 mG) in Na 7. This result justifies the conclusion of |7] that a 
weak resonance is preferable for formation of entangled atoms. The energy spectrum 
of non-condensate atoms 

{E, t) = mpusip, t)/{2n^), (19) 

presented in figure [^b; is rather narrow, just as in the case of Na [7]. 

In conclusion, this letter presents an application of the parametric approximation 
[7| to a description of BEC losses in the JILA ^^Rb experiments [Bj. The method, 
while disregarding inhomogeneities, offers a rather transparent demonstration of the 
relation of the excessive losses observed to the quantum many-body effects of Bose 
enhancement, and provides detailed information on the nature of the entangled atom 
pairs produced in the process. 
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